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Abstract. We analyze numerically the spin-dependent transport through 
coherent chains of three coupled quantum dots weakly connected to external 
magnetic leads. In particular, using the diagrammatic technique on the Keldysh 
contour, we calculate the conductance, shot noise and tunnel magnetoresistance 
(TMR) in the sequential and cotunneling regimes. We show that transport 
characteristics greatly depend on the strength of the intcrdot Coulomb 
correlations, which determines the spacial distribution of electron wave function 
in the chain. When the correlations are relatively strong, depending on the 
transport regime, we find both negative TMR as well as TMR enhanced above 
the Julliere value, accompanied with negative differential conductance (NDC) and 
super-Poissonian shot noise. This nontrivial behavior of tunnel magnetoresistance 
is associated with selection rules that govern tunneling processes and various 
high-spin states of the chain that are relevant for transport. For weak interdot 
correlations, on the other hand, the TMR is always positive and not larger than 
the Julliere TMR, although super-Poissonian shot noise and NDC can still be 
observed. 
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1. Introduction 

Tunnel magnetoresistance (TMR) is a measure of system transport properties change 
when the magnetic configuration of the device switches from parallel to antiparallel 
alignment pQ. The tunneling current is usually larger in the parallel configuration, 
when transport occurs between the majority-majority and minority-minority spin 
bands, than in the antiparallel configuration, where electrons tunnel between majority 
and minority spin bands, which gives rise to positive TMR effect. The TMR has 
been analyzed in various systems, including single-electron transistors and quantum 
dots 0[lll[3[a[l[Hl[a[I3[IIl[Il[Il[Il[IillIS]. Infact,a great deal of theoretical 
and experimental investigations has been devoted to spin-polarized transport through 
quantum dot structures. This is because quantum dots coupled to ferromagnetic leads 
are ideal candidates to study the fundamental interactions between spins and charges 
[12 [111 [19, 20 . Furthermore, such systems are also being considered for applications 
in future spintronic devices as well as for quantum computing [21] , However, most 
of existing theoretical considerations of spin-dependent transport in quantum dots 
involved only single and double dot systems [3ll[3[i[7l[i[S[T0l[Ill[Ta[T3[ll[ig[16], 
while experiments were carried out mainly for single dot structures [22J, [23], [24J [25J, [26] 
EB EH ESI [301 EH- In particular, it has been shown [7 that the TMR in quantum dots 
weakly coupled to ferromagnetic leads is generally smaller than the value given by the 
Julliere model pQ, TMR Ju11 = 2p 2 /(l - p 2 ), where p is the spin polarization of the 
leads, which is characteristic of tunneling through a single tunnel junction. This result 
is rather intuitive, as by embedding a quantum dot structure between ferromagnetic 
electrodes, the tunneling processes through the system become incoherent due to spin- 
flip processes and spin relaxation in the dot, leading to suppressed TMR. Because 
the magnitude of TMR is generally conditioned by the interplay of spin-dependent 
tunneling, spin accumulation and various spin states that mediate the current, one 
may expect that in the case of multi-dot structures, where some high-spin molecular 
states may form, the behavior of TMR will be modified as compared to that observed 
in the case of single and double dots. 

To prove the above statement, in this paper we address the problem of tunneling 
through chains of quantum dots, consisting of three coherent dots, weakly coupled to 
external ferromagnetic leads. Very recently transport properties of triple quantum dots 
have become a subject of intensive studies due to various interesting effects that emerge 
in such structures [321 [331 EH ES GH EH EH E3 HQl [SI S21 |43l Sll HQ SSI S3 SHI |49l [50] . 
In particular, triple dots enable the investigation of spin-entangled currents [32] . 
dark states [371 [H] , or various interference effects [331 ES] • Although nonmagnetic 
properties of multi-dot structures have already been addressed both theoretically and 
experimentally, very little is known about their magnetic transport properties [5 lj . 
The goal of this paper is therefore to discuss the spin-polarized transport through 
coherent triple quantum dots. In particular, by employing the real-time diagrammatic 
technique, we calculate the current, differential conductance, TMR and shot noise in 
both the sequential and cotunneling regimes. We show that transport characteristics 
strongly depend on the strength of the interdot correlations, which determines the 
spacial distribution of electron wave functions in the chain. In the case of strong 
Coulomb correlations, we find that the TMR may take values larger than the Julliere 
TMR, which is associated with tunneling through high-spin molecular states of the 
quantum dot system. Moreover, we also predict negative TMR, due to an increased 
tunneling current in the antiparallel configuration, associated with spin accumulation 
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in the chain. In addition, we show that these effects may be accompanied with super- 
Poissonian shot noise and negative differential conductance (NDC). On the other hand, 
in the case of weak interdot Coulomb correlations, the TMR is always positive and not 
larger than the Julliere TMR, while we still observe super-Poissonian shot noise in the 
Coulomb blockade regime and negative differential conductance. Although here we 
consider chains consisting of only three quantum dots, similar behavior may be also 
observed in longer chains where transport occurs through high-spin molecular states 
and is governed by various selection rules. 

The paper is organized as follows. In the second section we describe the 
Hamiltonian of the quantum dot chain and briefly discuss the method used in 
calculations. Section III is devoted to numerical results, where we first analyze the 
transport characteristics in the case of strong interdot Coulomb correlations and then 
proceed to discuss the transport behavior in the case of weak interdot correlations. 
Finally, the conclusions are given in Sec. IV. 

2. Theoretical description 

2.1. Model 

The schematic of chain consisting of three quantum dots coupled to ferromagnetic 
leads is shown in Fig. [T] It is assumed that the magnetizations of the leads are 
oriented collincarly, so that the system can be cither in the parallel or antiparallel 
magnetic configuration. The Hamiltonian of the system is given by 

H = -fflcad + Htun + ^chain , (1) 

where the first part corresponds to noninteracting itinerant electrons in the left (r = L) 
and right (r = R) lead, iJi cad = J2 r Sko- eVko-c^Crka-, where s r ka is the energy of an 
electron with the wave vector k and spin a in the lead r, and cj kcr (c^ko-) denotes the 
respective creation (annihilation) operator. The second term of Eq. ([!} accounts for 
the tunneling processes between the leads and the quantum dot chain, 

-fftun = ^2 ( iLC Lk<7 d l<7 + fe4ka rf 3<x + h.C.) , (2) 
kcr 

where t r denotes the tunnel matrix elements between the lead r and the respective dot 
and dj a destroys a spin-cr electron in the dot j (j = 1, 2, 3). Note, that the first dot is 
coupled to the left lead, while the third dot is connected to the right lead, see Fig. Q] 
The strength of the coupling of the quantum dot chain to the spin-majority (spin- 
minority) electron band of the rth lead is given by, r^H ' = 2ir\t r \ 2 pt^ = T r (l±p r ), 
where T r = (r++T~)/2, while pt ( ~ ] and p r are the spin-dependent density of states 
for majority (minority) spin band and spin polarization in the lead r, respectively. In 
the following we assume Tl = Tr = T/2 and p^ = pn = p. 

Finally, the last term of the Hamiltonian describes the chain consisting of three 
quantum dots which is given by 

#chain = e 3 n i° + U n J1 n 3l 

3 CT 3 

+ U'2_j ( n l<7 n 2a' + n2c,n z< j>) 
<j<j' 

+ 1 ( d l* d2 ° + d la d 3a + h.C.^j , (3) 
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Figure 1. (color online) The schematic of chain consisting of three single- level 
quantum dots connected to external ferromagnetic leads. The hopping between 
neighboring dots is denoted by t, U' and U are the inter-dot and intra-dot 
Coulomb correlation energies, while IY and denote the couplings to the left 
and right leads. The magnetizations of the leads are assumed to form either 
parallel or antiparallel magnetic configuration, as sketched in the figure. 

with riju = dj a dj a begin the particle number operator on dot j, while £j and Uj denote 
the single-particle energy and on-level Coulomb correlation in dot j, respectively. The 
third part of i? c i ia in corresponds to the inter-dot Coulomb interaction, whose strength 
is given by U', while t describes the inter-dot hopping. As we are interested in rather 
low bias voltage regime, it is justifiable to assume that the dot energy levels are 
independent of the bias voltage. This assumption has also been verified numerically 
and even if one assumes a voltage drop on the outer dots of the order of 10%, the 
current flowing through the system becomes only very slightly modified as compared 
that in the case of zero voltage drop. Therefore, for the sake of clarity of further 
discussion, we take the energy-independent dot levels. Furthermore, we also assume 
that the system is symmetric, i.e., Ej = e and Uj = U (j — 1, 2, 3). 




2.2. Method 

In order to calculate the spin-polarized transport through a chain of three coherent 
quantum dots in the sequential and cotunneling regimes, we employ the real-time 
diagrammatic technique [TJ I52[ 153] . It generally consists in a perturbative expansion 
of the density matrix of the system and the operators of interest (for example the 
current operator) with respect to the coupling strength T. Time evolution of the 
reduced density matrix is given by a sequence of irreducible self-energy blocks, , 
on the Keldysh contour, corresponding to various transition events between the many- 
body states |x) and |x') of the quantum dot chain. On the other hand, the full 
propagation of the reduced density matrix is given by the Dyson equation, which 
can be further transformed into a general kinetic equation. In a steady state the 
kinetic equation is simply given by (SP) X = T5 XXo , and enables the calculation of 
occupation probabilities P x for the system to be in a many-body state Here, P is 
the probability vector, while the matrix S is given by the self-energy matrix S with 
one arbitrary row xo replaced by (T, . . . ,T) due to normalization, ~^2 X P x = 1. The 
current flowing through the system can be then found from [52 

7=-^Tr{S I P}, (4) 

where S 1 denotes the modified self-energy matrix S so as to take into account the 
number of electrons transferred through the system. 

To calculate the transport properties order by order in tunneling processes, we 
expand the self-energy matrices, £ = £ (1 ) + £ (2 ) + . . ., S 1 = S^^+S 1 ^^. . ., and the 
occupations, P = p(°) + P^ 1 ) + . . ., respectively. The self-energies in respective order 
can be then calculated using the corresponding diagrammatic rules [52] . The first 
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order of expansion corresponds to the sequential tunneling, whereas the second one to 
cotunneling. In this analysis we have calculated the self-energies up to the second order 
of perturbation series, so that we are able to resolve transport properties both in the 
sequential and cotunneling regimes [15] . The sequential tunneling dominates transport 
above a threshold voltage and is exponentially suppressed in the Coulomb blockade 
regime |54j . In the blockade regime, on the other hand, the dominant contribution to 
the current comes from cotunneling processes |55) . which take place through virtual 
states of the system and are only algebraically suppressed in the Coulomb blockade. 
As the influence of cotunneling on transport for bias voltages above the threshold for 
sequential tunneling is rather minor, the inclusion of second-order processes is crucial 
for a proper description of transport in the blockade regimes. 

In addition, in the following we will also analyze the zero-frequency current noise 
[56], S = dt((/(i)i(0) + I(0)I(t)) - 2(/) 2 ), where I is the current operator, 

I = (Ir - 4)/2, with J L(R) = -i(e/h)t L{R) Ek<T(4(R)k<r d i(3) ff - 4(3) CT c L(R)ka) being 
the current flowing from the first (third) dot to the left (right) lead. The formula for 
current noise derived within the real-time diagrammatic technique can be found in 
Ref. [53]. 

3. Numerical results 

In the following we will discuss the numerical results on the current, differential 
conductance, tunnel magnetoresistance and the shot noise of a chain of tree coherent 
single-level quantum dots in both the linear and nonlinear response regimes. Transport 
characteristics of such systems strongly depend on the internal parameters, in 
particular, on the ratio between inter-dot Coulomb repulsion U' and the hopping 
between the dots t, provided that U > U', \t\. The ratio can be tuned experimentally 
for example by changing the height of the barrier between the dots [57]. When the 
inter-dot Coulomb correlations are relatively strong, U'/\t\ > 1, the electrons in the 
ground state of the chain will be mostly occupying the outermost dots. On the other 
hand, for weak inter-dot Coulomb interactions, U'/\t\ < 1, this tendency will not 
be observed. Thus, depending on U'/\t\, the spacial distribution of the many-body 
chain states may become strongly modified. In this paper we will therefore discuss the 
transport characteristics in the two above mentioned situations. Furthermore, we also 
note that due to many intrinsic parameters of the system, there is a variety of transport 
regimes where different transport behavior can be observed. In the following, we will 
thus present general density plots, however only most interesting transport features 
will be discussed in greater detail. 

3.1. The case of strong inter-dot Coulomb interactions 

Figure [2] shows various transport characteristics of a quantum dot chain as a function 
of bias voltage V and the position of the dots' levels e. Because the position of the 
levels can be experimentally changed by sweeping the gate voltage, Fig. [2] effectively 
presents a bias and gate voltage dependence of transport characteristics. The total 
(first plus second order) differential conductance for the parallel (Gp) and antiparallel 
(Gap) magnetic configurations is shown in Fig. [2]Ja) and (b), respectively. First of all, 
one can see that the differential conductance displays characteristic Coulomb diamond 
pattern, with Coulomb blockade regimes at low transport voltages. By lowering the 
position of the dots' levels the chain becomes consecutively occupied with electrons. 




Figure 2. (color online) The total (first plus second order) differential 
conductance in the parallel (a) and antiparallel (b) alignment, the total TMR 
(c), the TMR calculated using only first-order processes (d) and the total Fano 
factor for parallel (e) and antiparallel (f) configuration of the system. The figures 
were calculated for the case where the inter-dot Coulomb interaction is larger than 
the inter-dot hopping, U' /\t\ > 1. The parameters are: U 1 = 4 meV, t = —2 meV, 
U = 10 meV, k B T = 0.15 meV, T = 0.1 meV and p = 0.5. The Fano factor at 
low bias voltages diverges, therefore this transport regime is marked with white 
stripe in (e) and (f). 



In the case considered here, the quantum dot chain can accommodate up to six 
electrons, i.e. each dot can be doubly occupied. For such values of e when the two 
neighboring charge states become degenerate, there is a peak in the linear conductance. 
On the other hand, with increasing the bias voltage, out of the Coulomb blockade 
regime, there are additional lines visible in the differential conductance associated with 
tunneling through excited states of the system. Furthermore, due to the contact with 
ferromagnetic leads, the tunneling processes in the system become spin-dependent 
and, consequently, transport depends on the magnetic configuration of the system. 
In the parallel configuration the majority (minority) electrons of one lead tunnel 
to the majority (minority) spin band of the other lead, whereas in the antiparallel 
configuration they tunnel to the minority (majority) spin band. This is why the 
conductance in the antiparallel configuration is generally suppressed as compared to 
that in the parallel configuration, Gp > Gap, see Fig. [^a) and (b). This difference 
in turn gives rise to nonzero tunnel magnetoresistance which is plotted in Fig. [He) . 

The TMR reflects the change of system transport properties when switching the 
magnetic configuration of the device from parallel to antiparallel one. It is qualitatively 
defined as [D [2 S], TMR = I P /I A p - 1, where I P (7 A p) it the current flowing 
through the system in the parallel (antiparallel) magnetic configuration. Usually 
the conductance in the parallel configuration is larger than that in the antiparallel 
one, giving rise to positive TMR. In particular, for a single ferromagnetic tunnel 
junction the TMR can be described by the Julliere model lj, TMR Ju11 = 2p 2 /(l - p 2 ) 
(TMR Ju11 = 2/3 for p = 0.5 assumed in calculations). Intuitively, one may expect that 
placing a quantum dot molecule between the two leads (where tunneling processes are 
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generally incoherent) will decrease the TMR. This is in fact what is observed in most 
quantum dot structures - for symmetric systems and in the absence of magnetic field, 
the TMR in the weak coupling regime is positive and not larger than TMR Ju11 pj [TB] . 
In the case of tunneling through quantum dot chains considered here, however, we 
predict a nontrivial behavior of the TMR. Depending on the transport regime, we 
find both the TMR enhanced above the Julliere value as well as negative TMR effect, 
see Fig. [2jc) . The mechanisms responsible for these effects will be discussed in more 
detail in the sequel. 

We also note that the TMR is directly related to the ratio between the currents in 
the two magnetic configurations, so that its magnitude does not necessarily dependent 
on the magnitude of tunneling current. This makes the TMR a vary sensitive quantity 
for analyzing transport properties in various regimes, especially where sequential 
tunneling is suppressed due to Coulomb correlations and transport occurs mainly 
through higher-order tunneling processes. For comparison, in Fig. HJd) we have also 
plotted the TMR calculated using only the first-order tunneling processes, TMR' 1 '. 
It can be seen that cotunneling processes modify the TMR mainly in the blockade 
regimes and give rise to strong dependence of TMR on the occupation number of the 
quantum dot chain. On the other hand, out of the blockade regimes, the current is 
predominantly mediated by sequential tunneling and one finds that TMR and TMR* 1 ' 
become comparable, although not equal. 

In addition, we have also calculated the Fano factor, F — S/S p , in both magnetic 
configurations, see Fig. [2{e) and (f). The Fano factor describes the deviation of the 
shot noise from its Poissonian value, S p = 2e\I\, which is characteristic of uncorrelated 
tunneling. When transport is mediated only by elastic cotunneling processes, the noise 
is Poissonian, F — > 1, however, once the spin-flip cotunneling is allowed the noise 
can be enhanced to become super- Poissonian, F > 1, due to bunching of inelastic 
processes [58 . Furthermore, in the sequential tunneling regime, transport is mainly 
dominated by Coulomb correlations which decrease the noise and the Fano factor 
is generally sub-Poissonian, F < 1 59J. This behavior can be in fact observed in 
Fig. [He) and (f), where in the cotunneling regime the Fano factor can take large 
super-Poissonian values, while in the sequential tunneling regime it becomes rather 
suppressed. It can be also seen that the general behavior of the Fano factor in the 
parallel (Fp) and antiparallel (Fap) magnetic configurations is quite similar, although 
the magnitude of the noise is larger in the antiparallel configuration. On the other 
hand, in the low bias voltage regime, the noise is dominated by thermal noise while 
the current tends to zero, which leads to a divergency in the Fano factor. Therefore 
this transport regime is marked with white stripes in Fig. \2le) and (f). 

3.1.1. Linear response regime In the linear response regime, the transport behavior is 
mainly conditioned by the ground state of the system and its evolution when changing 
the position of the dots levels. The ground state energies Eq s together with respective 
quantum numbers of states {Q, S} are shown in Table [T] with 

Q = n i' J ' S = 2 d\ a 3 aa id ja i , (5) 

ja jaa' 

denoting the total charge and total spin of the quantum dot chain. Because in the 
absence of external magnetic field the Hamiltonian of the chain, Eq. (|3|), commutes 
with Q and S 2 , one can solve the eigenvalue problem by diagonalizing H chain in 
respective blocks {Q,S}. Furthermore, by using the full spin SU(2) symmetry the 
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Figure 3. (color online) The linear conductance (a) in the parallel (solid line) and 
antiparallel (dashed line) magnetic configurations and linear response TMR (b) as 
a function of the dots' level position s. The parameters are the same as in Fig. [2] 
The dotted line in (b) shows the linear TMR calculated by using only first-order 
processes - sequential TMR is constant and given by TMR' 1 ' = p 2 /(l — p 2 ). 



size of the Hilbert space is effectively reduced from 64 to 35 multiplets, which may 
be crucial for analytical discussion of decoupled quantum dot chain. In numerical 
calculations, however, we have used the 64-state space with the following many-body 
states: \x) = 1X1X2X3); where Xj — 0, 1\ |, d denotes the state with zero electrons, one 
spin- up, spin-down electron and two electrons on the dot j. This is because in the 
case of the spin-dependent coupling to ferromagnetic leads, the Hamiltonian of the 
whole system possesses only the S z symmetry. In Table Q] we show the corresponding 
quantum numbers, the dimension of the Hamiltonian blocks, degeneracy of states and 
ground state energies of the decoupled quantum dot chain. It turns out that the 
ground state energies Eg s for smaller blocks can be easily calculated, however, for 
larger blocks the formulas become too lengthy to be presented here. We thus list the 
explicit expressions for low-dimension Hamiltonian blocks, while for the other blocks 
we just state which energy is the lowest one in respective charge sector Q. The explicit 
matrices for Hamiltonian blocks Hq^$ together with the definition of states for total 
charge and total spin symmetries can be found in the Appendix. 

The linear conductance as well as the total TMR are shown in Fig. 02a) and 
(b). The linear conductance displays characteristic resonance peaks whenever two 
neighboring charge states become degenerate. The resonance energies can be estimated 
from Table[T]by solving min{i?Q +1 s } — min{i?Q s } = 0, where one needs to take the 
minimum energy for given Q. The conductance in the parallel configuration is larger 
than the conductance in the antiparallel configuration, see Fig. E2a), which results in 
positive linear TMR, see Fig. EHb). For comparison we have also plotted the TMR 
obtained using only first-order tunneling processes, which is constant and given by 
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Table 1. The charge Q and the total spin S of three quantum dots coupled 
in series, dimension of respective {(3,5} blocks, D, degeneracy of states, d, and 
ground state energies, Sq S - The analytical formulas for ground state energies 
of larger blocks are rather lengthy, therefore here we only state which states have 
lower energy for given Q. The explicit matrices for Hamiltonian blocks are given 
in the Appendix. 
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TMR* 1 ' = p 2 /(l -p 2 ), see the dotted line in Fig.^b). The total TMR, on the other 
hand, shows a nontrivial dependence on the position of the dots' levels e. As shown in 
the case of single quantum dots [3 [16] , the magnitude of linear TMR is directly related 
to the type of cotunneling processes that drive the current in respective transport 
regimes. Among various cotunneling events, one can distinguish processes that affect 
the magnetic state of the quantum dot system (inelastic spin-flip processes) and the 
ones that do not affect the quantum dot chain (elastic non-spin-flip processes). In 
the case when each dot of the chain is either empty or doubly occupied, only elastic 
processes are possible, however, in the other cases, the spin-flip processes become also 
allowed. 

When the chain is empty (Q = 0) or fully occupied (Q = 6), TMR is maximal and 
equal to the Julliere value, TMR = 2p 2 /(l - p 2 ), see Fig. [3(b). This is related to the 
fact that in these transport regimes only the non-spin-flip cotunneling processes are 
allowed and transport is fully coherent - the co-tunneling electrons are not scattered 
at the chain at all. Interestingly, also for Q = 4 the TMR becomes equal to TMR , 
which indicates that only elastic cotunneling contributes to the linear conductance. In 
fact, for Q = 4 the ground state of the chain is S — 0, see Table [TJ Furthermore, for 
assumed parameters, i.e. in the case of strong interdot correlations, U'/t > 1, it turns 
out that in the ground state the chain is occupied with two electrons in the outermost 
dots, so that the ground state is |d0d). In this case only the non-spin- flip cotunneling 
is allowed, which yields the maximum TMR. However, the situation changes once the 
spin-flip processes become possible, which happens in the other transport regimes. In 
particular, for Q = 1, 2, 3, the linear TMR becomes suppressed to approximately a half 
of TMR Ju11 and its dependence on e is rather complex. In these regimes the current 
is mainly mediated by inelastic spin- flip cotunneling. For Q — 5, on the other hand, 
the TMR becomes slightly enhanced, although it is still lower than TMR Ju11 . Because 
for Q = 5 the ground state is a doublet |dcxd), transport is due to both elastic and 
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Figure 4. (color online) The bias voltage dependence of the current (a), 
differential conductance (b), Fano factor (c) in the parallel (solid line) and 
antiparallel (dashed line) magnetic configurations and the resulting TMR (d) for 
e = —7.5 meV. The other parameters are the same as in Fig. [2] and Iq denotes 
the maximum current given by Iq = eT/h. 



inelastic processes - the former (latter) ones tend to increase (decrease) the TMR, 
so that the magnitude of TMR is between the values corresponding to Q — 1,2,3 
and Q — 0,4, 6 transport regimes. Finally, we also note that at resonances the total 
TMR drops to the value approximately given by TMR (1) , as for resonant energies the 
first-order processes become possible and are dominant. 
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3.1.2. Enhanced TMR and negative differential conductance The bias dependence of 
the current, differential conductance, Fano factor and TMR is shown in Fig. [U The 
transport characteristics were calculated for s = —7.5 meV, see also Fig. [51 which 
corresponds to the case when at equilibrium the quantum dot chain is in the spin 
doublet state, {Q = 3, S — 5}, so that the ground state is doubly degenerate. With 
increasing the bias voltage, more and more states start participating in transport 
and the current increases. However, it can be seen that the bias dependence of the 
current is not monotonic - after the first Coulomb step, the current starts to decrease 
with V, leading to negative differential conductance (NDC), which is present in both 
magnetic configurations, see Fig. BJb). The suppression of the current is associated 
with selection rules that govern the respective sequential transitions, i.e. only the 
transitions that change the total charge of the chain by 1 and the total spin by \ 
are allowed. When raising the transport voltage, the following excited states {Q, S} 
become active in transport: {2,0} singlet and {2, 1} triplet, then {3, \} doublet, and 
then {3, |} quadruplet, respectively. Although the sequential transitions between 
Q = 3 doublet and Q — 2 singlet and triplet are possible, the transitions between 
Q = 3 doublet and Q — 3 quadruplet are prohibited as they obey neither the charge nor 
the spin selection rules. In fact, once the system gets trapped in the quadruplet state, 
transport becomes suppressed, see Fig. \M,&) for V ~ 5 mV, which leads to negative 
differential conductance. This is because transitions involving {3, |} state can occur 
only through tunneling from/into the triplet state {2, 1}. In the case of U/\t'\ > 1, 
the two electrons in the triplet state are localized in the outermost dots, while in the 
quadruplet state the three electrons are distributed uniformly between the dots. Thus, 
for a transition between {3, §} and {2, 1} to occur, one needs to put or remove an 
electron from the middle dot, which is however suppressed as (t/U') 2 . On the other 
hand, tunneling processes in which the state of the system changes between the Q = 3 
doublet and Q — 2 singlet are rather independent of the ratio t/U', as they can occur 
through the outmost dots. Consequently, transitions involving the quadruplet state 
are relatively slow, while the other ones are much faster. The competition between 
such slow and fast transport channels may in turn lead to large current fluctuations 
[3"5] . This can be seen in Fig. 0Jc) , where for voltages corresponding to the transport 
regime where the current is suppressed, super-Poissonian shot noise is observed. When 
increasing the bias voltage further, more excited states become available for transport 
and the current starts increasing again, see Fig. Ufa), while the shot noise becomes 
suppressed to sub-Poissonian value, which is typical of charge-correlated sequential 
transport, see Fig. H^c). 

In the antiparallel configuration, on the other hand, due to the asymmetry of 
tunneling processes between the left and right leads, there is a nonequilibrium spin 
accumulation in the chain. For positive bias voltages, the occupation probability of 
highest-weight spin states is much increased as compared to the other components 
of particular state. This is because the spin-up electrons tunneling from the left 
lead to the chain and the spin-down electrons tunneling out of the chain to the 
right lead belong to the majority-spin bands, and the positive spin component 
becomes accumulated in the chain. Consequently, fewer states are available for 
transport as compared to the parallel configuration, so that in the antiparallel 
configuration the transport channel involving the quadruplet state {3, |} becomes 
even less transmitting. This leads to several interesting features. First of all, the 
current is more suppressed in the antiparallel configuration than in the parallel one 
which leads to an enhanced NDC, see Fig.[4|b). O n the other hand, this more effective 
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current suppression also reflects itself in an enhanced super-Poissonian shot noise, 
see Fig. HJc). Furthermore, for voltages where the current is suppressed we observe 
the TMR enhanced above the Julliere value. The TMR in this transport regime 
is approximately given by TMR w |TMR J , see Fig. Sfd). As shown in previous 
considerations, such enhancement of TMR in serial quantum dots weakly coupled to 
external leads can occur mainly in asymmetric systems or in the presence of external 
magnetic field [HI [16]. Here, we observe TMR > TMR Ju11 in the absence of magnetic 
field and for fully symmetric system. This is just a generic feature of transport through 
chains of quantum dots, where due to selection rules the system may be trapped in 
some high-spin states. 

3.1.3. Negative TMR and super-Poissonian shot noise Another interesting transport 
behavior can be observed in the case where the chain is in the ground state with four 
electrons. The current, differential conductance, shot noise and TMR as a function of 
the bias voltage for e = —15 meV are shown in Fig. [5] Due to large interdot Coulomb 
correlations, U'/\t\ > 1, the ground state is non-degenerate, with doubly occupied 
outermost dots, | dOd) . The nearest excited states are respectively: doublets {3,^} 
and {5, i}, and triplet {4, 1}. These states are relatively close to each other and start 
taking part in transport for voltages around the threshold for sequential tunneling. 
At low bias voltages the system is in the Coulomb blockade, see Fig.[5{a) and (b), and 
transport is due to elastic cotunncling processes, which yield the Poissonian shot noise 
and TMR given by the Julliere value, see Fig.[5]Jc) and (d). The situation changes once 
the transport voltage approaches the threshold, V f=s 5 mV, then the TMR suddenly 
drops and changes sign, while the shot noise becomes strongly enhanced. This is 
associated with tunneling processes that become allowed in this transport regime. 
First of all, the inelastic cotunneling processes become possible for each doublet state, 
i.e. {3, j} an d {5, \}. Furthermore, around the threshold voltage, the sequential 
processes also start participating in transport. The first-order transitions occur first 
between the Q = 4 singlet and Q — 3, Q — 5 doublets. It is worth noting that the 
spacial distribution of the wave function is different for these two doublets. For Q = 3 
the electrons are equally distributed between the three dots, while for Q = 5 the 
outermost dots are fully occupied while the middle dot is singly occupied to minimize 
the interdot correlations. Consequently, the transport channel involving the state 
{5, i} is slower than that involving the state {3, i}, similarly as in the case discussed 
in previous subsection. It turns out that the interplay of various first and second-order 
tunneling processes, where particular events occur at different rates, which exists for 
transport voltages around the threshold for sequential tunneling, leads to large current 
fluctuations. As a result, we observe an enhanced super-Poissonian shot noise in 
both magnetic configurations of the system, see Fig. OUc). On the other hand, when 
the voltage increases further, sequential processes dominate transport and the noise 
becomes generally sub-Poissonian. 

An interesting transport feature visible around the threshold voltage is the 
negative TMR effect, see Fig. [5jd). To understand this behavior, one needs to 
realize a very delicate difference between probability distributions in the two magnetic 
configurations. In the Coulomb blockade regime the chain is in the singlet state |dOd) 
with probability equal to unity, irrespective of magnetic configuration of the system. 
However, once the bias voltage approaches the threshold voltage, the occupation 
probability of states {3, |}, {5, 4}, and {4, 1} starts slowly increasing. In addition, 
it turns out that in the antiparallel configuration the probability of the highest- 
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Figure 5. (color online) The bias voltage dependence of the current (a), 
differential conductance (b), Fano factor (c) in the parallel (solid line) and 
antiparallel (dashed line) magnetic configurations and the resulting TMR (d) for 
e = — 15 meV and the parameters as in Fig. [2] 

weight quadruplet state is also nonzero, and it is slightly larger than the occupation 
probabilities of the above-mentioned doublets and triplet. The enhanced occupation 
probability of {3, |} is purely due to nonequilibrium spin accumulation, it is thus not 
present in the parallel configuration. This is in fact what is crucial for the occurrence 
of negative TMR. In the antiparallel configuration the current can in addition flow due 
to cotunneling and thermally-activated first-order transitions involving the quadruplet 
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Figure 6. (color online) The total differential conductance in the parallel (a) and 
antiparallel (b) alignment, the total TMR (c), the TMR calculated using only first- 
order processes (d) and the total Fano factor for parallel (e) and antiparallel (f) 
configuration of the system for the case when the inter-dot Coulomb interaction 
is smaller than the inter-dot hopping, U' /\t\ < 1. The parameters are: V = 0.5 
meV, t = -2 meV, U = 10 meV, k B T = 0.15 meV, T = 0.1 mcV and p = 0.5. 
The Fano factor at low bias voltage diverges, therefore this transport regime is 
marked with white stripe in (e) and (f). 



state, which is not possible in the parallel configuration. As a result, the current in the 
antiparallel configuration becomes larger than the current in the parallel configuration, 
yielding negative TMR effect. This can be also seen in the differential conductance 
where the first peak in the antiparallel configuration occurs at slightly lower bias 
voltage than in the parallel one, see Fig. [5th). With increasing the bias voltage further, 
the excited states start participating in transport and the system apparently exhibits 
a normal spin- valve behavior [6l [7], with the current in the parallel configuration 
larger than in the antiparallel one and, thus, with positive TMR, see Fig. [Sfd). 
It is interesting to note that the increased occupation probability of the highest- 
spin component of the quadruplet state in the antiparallel configuration was also 
responsible for the enhanced TMR effect discussed in previous subsection, whereas 
here it lead to negative TMR. The negative TMR however occurs on the edge of 
the Coulomb blockade regime, while the enhanced TMR develops in the sequential 
tunneling regime. 

3.2. The case of weak inter-dot Coulomb interaction 

The differential conductance, TMR and Fano factor in the case of U'/\t\ < 1 are 
shown in Fig. [6l For U'/\t\ < 1, the electrons in particular states are distributed 
rather uniformly over the three dots, contrary to the previous case where electrons 
where localized in the outermost dots to minimize the Coulomb correlation energy. 
This results in a more symmetric behavior of transport characteristics with respect to 
the middle of the Coulomb blockade regime with Q — 3, which is due to particle- hole 
symmetry. Moreover, as most of the effects observed U'/\t\ > 1 were mainly associated 
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with spacial distribution of the wave function, one may expect their strong dependence 
on the ratio U'/\t\. This is in fact what can be observed. For example, when U'/\t\ < 1, 
the negative TMR and TMR enhanced above the Julliere value are not present, 
although super-Poissonian shot noise and negative differential conductance can still 
be found in some transport regimes. Of course, the difference between transport 
characteristics in the case of weak and strong interdot correlations reveals itself mainly 
in transport regimes where the states with more than a single electron become relevant, 
see Figs. [2] and El 

The differential conductance in the parallel and antiparallel configurations is 
shown in Fig.[6{a) and (b), respectively. First of all, we note that because the energy of 
interdot correlations is now changed, the Coulomb diamond structure is different from 
that shown in Fig. [2] The largest Coulomb diamond develops for Q — 3, while the 
size of other diamonds is much decreased, except for empty and fully occupied chain, 
see Fig. HJa) and (b). The conductance in the parallel configuration is larger than in 
the antiparallel one and the TMR is positive in the whole range of bias voltage V and 
the level position e, see Fig.[6tc). For comparison, the TMR calculated using only the 
sequential tunneling processes is shown in Fig.[6]^d). The main difference between the 
density plots for the total TMR and TMR (1) can be seen in the Coulomb blockade 
regimes where cotunneling dominates the current. It can be seen that the total TMR in 
the linear response regime displays nontrivial dependence on the occupation number of 
the chain. For empty and fully occupied chain, TMR is given by the Julliere value and 
it is much suppressed in other blockade regimes due to spin-flip cotunneling processes. 

In the nonlinear response regime when at equilibrium the chain was in the charge 
state Q — 3, there is reminiscent of effects found in the case of U'/\t\ > 1, see Fig. |U 
Now, one can also observe an enhanced TMR, although its magnitude is slightly lower 
than the Julliere TMR. This enhanced TMR is accompanied with negative differential 
conductance and super-Poissonian shot noise, which are more visible in the antiparallel 
configuration. The mechanism leading to these effects is similar to that discussed in 
the previous subsection and is mainly associated with transport channel involving the 
quadruplet state and spin accumulation in the antiparallel configuration. 

Furthermore, in the nonlinear response regime of the Coulomb blockade regime 
with Q = 2, the TMR becomes suppressed, being very close to zero. In this transport 
regime the chain at equilibrium is in the singlet state {2,0}, and the excited states 
are consecutively {2,1} and {1,^}- With increasing the bias voltage, occupation 
probability of excited states starts increasing. It turns out that in the parallel 
configuration all the components of Q = 1 doublet and Q = 2 triplet are relevant, 
while in the antiparallel configuration, due to spin accumulation, only the highest- 
weight components, however with slightly larger occupation probabilities. This leads 
to an increased current in the antiparallel configuration, so that the currents in the 
two configurations become roughly comparable, yielding very small TMR. In addition, 
with raising the bias voltage the shot noise becomes enhanced and reaches maximum 
for voltages around the threshold for sequential tunneling, which is associated with 
bunching of inelastic cotunneling processes. Similar behavior can be also observed 
for the Coulomb blockade regime with four electrons in the chain, which is due to 
particle-hole symmetry. It is also interesting to note that the behavior of transport 
characteristics in the cotunneling regime with Q = 2 is very weakly affected by the 
ratio U'/\t\, see Figs. [5] and |B] However, for Coulomb blockade regimes with more 
electrons in the ground state, transport properties become completely modified due 
to different spacial distribution of wave functions, compare for example the Coulomb 
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blockade regime with Q = 4 in Figs. and [51 

Finally, we also note that in the case of weak interdot correlations the shot 
noise is rather sub-Poissonian in the whole sequential tunneling regime, irrespective of 
magnetic configuration of the system. On the other hand, the super-Poissonian shot 
noise is only found in the Coulomb blockade regimes where the chain is the charge state 
with two, three or four electrons, which is due to bunching of inelastic cotunneling 
processes. 

4. Conclusions 

Wc have analyzed the linear and nonlinear transport properties of chains of quantum 
dots consisting of three coherent single-level quantum dots weakly coupled to 
external ferromagnetic leads. By employing the real-time diagrammatic technique, 
we have calculated the current, differential conductance, shot noise and tunnel 
magnetoresistance in the case of strong (U' /\t\ > 1) and weak (?7'/|t| < 1) interdot 
correlations. By changing the ratio U'/\t\, one can effectively change the spacial 
distribution of electron wave functions of the chain. When U'/\t\ > 1, the electrons 
tend to be localized in the outermost dots, while for U'/\t\ < 1 the electrons are 
distributed rather uniformly over the dots. 

In particular, in the case of large interdot correlations we have shown that 
the TMR strongly depends on the transport regime and can take negative values 
as well as values exceeding the TMR given by the Julliere model. The enhanced 
TMR occurs in the nonlinear response regime when the chain is occupied by three 
electrons at equilibrium and is associated with a suppressed current in the antiparallel 
configuration due to trapping of the quantum dot chain in some high-spin states. In 
addition, the suppression of the current gives rise to negative differential conductance 
and super-Poissonian shot noise. We have also shown that the TMR may change sign 
and become negative. This happens in the cotunneling regime with four electrons in 
the chain when the bias voltage approaches threshold voltage for sequential tunneling. 
The negative TMR is then associated with increased tunneling through the highest- 
weight spin state (quadruplet) of the chain. Furthermore, we have also shown that 
negative TMR is accompanied with large super-Poissonian shot noise due to the 
interplay between various inelastic cotunneling and sequential processes that start 
contributing to the current around the threshold voltage. On the other hand, when 
the interdot correlations are weak, most of the effects found in the case of U' /\t\ > 1 
become smeared out. In particular, the negative TMR and TMR enhanced above the 
Julliere value are not present, although super-Poissonian shot noise and the negative 
differential conductance can still be observed. 

Finally, we note that although the results presented here were calculated for 
chains of three quantum dots, similar behavior may be in principle observed for longer 
chains, where transport is governed by various selection rules and the current can flow 
due to tunneling through high-spin molecular states of the system. 
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Appendix 

The localized basis of quantum dot chain states is defined as \x) — IX1X2X3), 
where Xj = 0, j,J.,d denotes zero electrons, one spin- up, spin-down electron and 
doubly occupied dot j. Using the SU(2) symmetry for total spin one can reduce 
the Hilbcrt space from 4 3 states to 35 multiplcts. The reduction of Hilbert space 
is important rather for analytical calculations, while in numerical calculations we 
have used the basis of 64 states. In the following, we give the explicit matrices 
for blocks of the Hamiltonian i/ c hain m the basis of total charge Q and total spin 
5, where to define the spin SU(2) basis we have taken the highest-weight spin 
states. In the block {Q = 2, S = 0}, the states are: Sf° = |d00), S^° = |0d0), 
S 2 3 >° = |00d), = ^(| U 0) - I IT 0)), 5 2 2 3 ° = ^(|0 U) - |0 IT)), and 

Sis ~ "^f(l T I) — | | t)), respectively. The Hamiltonian block in this basis 
is given by 



/ 


£21 








V2t 








\ 







£21 





V2t 


V2t 













£21 





V2t 









V2t 


V2t 





£22 





t 









V2t 


V2i 





£22 


t 




V 











t 


t 


£23 


/ 



(•1) 



where £21 = 2e + U, £22 = 2e + U', and e 2 3 = 2s. The block of i? c hain for 
{Q = 4, S = 0}, Hifi, has similar structure to #2,0 due to particle hole symmetry, 
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where £41 = 4e + 2U + 4U', £ 42 = 4£ + U + 3U', and e 43 = 4e + U + W. The 
states in block {Q = 4, S — 0} are explicitly given by: S*^ = |0dd), S^' = |d0d), 
S$>° = |dd0), S 4 12 ° = ^(| Tl d) - I |T d)), S&° = ^(|d Tl) - |d IT)), and 
Sis = "^7f(ITdJ.) — j j d t)), respectively. On the other hand, the Hamiltonian 
block for {Q = 3, S — |} is the largest one, with the states defined as follows: 
= |d T 0), = |0 T d), = I t d0), = |0d |), £)^' = | | Od), 

= |do t), = ^(| TIT) - 1 ITT)), and D 3 /* = J\\ TTI) - ^(| TIT) + 1 ITT)), 
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respectively. The block H 3 1 is explicitly given by 
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3e + U, and e 33 = 3e + 2U'. For completeness, we 



where £31 = 3e + U + 2U' ', £32 

also give the matrices for smaller blocks of the Hamiltonian, 

/ £ t 

\ * e 

/ 2e + U' 
H2.1 = t 

V 

The states for block {Q = 1,S= i} are: d\'* = \ | 00), ' = |0 | 0), D3' = |00 |), 




(■4) 



(.5) 



TT 0), T 2 ia = I T T>, 

and ^4.1 have similar 



while for block {Q = 2.5' = 1} they are given by: T l ' = 
= |0 TT); respectively. The Hamiltonian blocks H 5 
structure to blocks H x 1 and #2,1 due to particle hole symmetry. On the other hand, 
blocks -ffo.O: H 3 3 and ^6,0 are trivially one dimensional, see also Table [TJ 
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